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ABSTRACT 

We present the orbital and physical parameters of the detached eclipsing binary VI200 Cen¬ 
tauri (ASAS J135218-3837.3) from the analysis of spectroscopic observations and light curves 
from the All Sky Automated Survey (ASAS) and SuperWASP database. The radial velocities 
were computed from the high-resolution spectra obtained with the OUC 50-cm telescope and 
PUCHEROS spectrograph and with 1.2m Euler telescope and CORALIE spectrograph us¬ 
ing the cross-correlation technique TODCOR. We found that the absolute parameters of the 
system are Mi = 1.394 ± 0.030 M 0 , M 2 = 0.866 ± 0.015 Mq, Ri — 1.39 ± 0.15 R©, 
i?2 = 1-10 ± 0.25 Rq. We investigated the evolutionary status and kinematics of the binary 
and our results indicate that VI200 Centauri is likely a member of the Hyades moving group, 
but the largely inflated secondary’s radius may suggest that the system may be even younger, 
around 30 Myr. We also found that the eclipsing pair is orbited by another, stellar-mass object 
on a 351-day orbit, which is unusually short for hierarchical triples. This makes V1200 Cen a 
potentially interesting target for testing the formation models of multiple stars. 

Key words: binaries; eclipsing - binaries: spectroscopic - binaries: pre-main sequence - stars: 
fundamental parameters - stars; individual; VI200 Cen 


1 INTRODUCTION 

Detached, double-lined spectroscopic binaries that are also eclips¬ 
ing provide an accurate determination of stellar mass, radius and 
temperature for each of their individual components, and hence 
constitute a strong test of sing le star stellar evolution theory 
jUastennet & Vails Gabaudl2002n . Eclipsing binary (EB) stars are 
very important for the study of stellar astrophysics. Their partic¬ 
ular geometrical layout, their dynamics and radiative physics en¬ 
able a detailed and accurate modelling and analysis of the acquired 
data, and all ow to measure many ba sic physical parameters of the 
components iHelminiak et alj|2009h . Studies of binary stars by all 
the techniques available in modem astrophysics allow to measure 
a wide range of parameters for each of the component stars, with 


* Based in part on observations collected through through CNTAC propos¬ 
als CN-2012A-21, CN-2012N-36 and CN-20I3A-93. 
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some of them determined with very high accuracy (e.g., uncertain¬ 
ties of less than 1%). 


This also refers to the pre-main-sequence (PMS) stars. 
Through a complete analysis of spectroscopy and photome¬ 
try of these systems orbital and physical parameters of the 
two stars can be accurately derived. However, until recently 
there were only seven known low-mass pre- main sequence EBs 
with M < 1.5 Mp,: RXJ 0529.4+ 0041A i lCovino et al.l l2000l. 
2004h . V117 4 Ori jStassun et al.l2004. 2 MJ0535-05 lls tassun_et_^ 


2006[ |2007h. JW380 jirwin et alj|2007h . P ar 1802 llCargile et al 


2008l:[Stassun et al 


and MML53 jHebb et al.ll2010l) . 


ASAS J0528-H03 dStempels et alj2008h . 


Mutiple systems of two or more bodies seem to be the norm 
at all stellar evolutionary stages, according to observations, and 
it has been accepted for some time now that binarity and mul- 
tiplicity is establishe d as the principal channel of star formation 
jReipurth et alJl20T4ll . Eor this reason, knowledge of the formation 
of multipl e stars becomes necessary to understand star formation 
in general jPelgado-Donate ^ aj2004ll . Moreover, observations of 
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young multiple systems allow to test evolutionary models in early 
stages of stellar evolution. 

On the other hand, it is important to increase the sample of 
young stars (especially PMS) because main-sequence, solar-type 
stars are well described by stellar evolution models (observations 
agree well with theoretical isochrones), but recent measurements of 
the stellar properties of low-mass dwarfs and young PMS stars re¬ 
main problematic for the existing models dMorales-Calderon et al.l 

l2012h . 


In this paper we present the orbital and physical analysis 
of an eclipsing binary V1200 Centauri (HD 120778, HIP 67712, 
ASAS JI352I8-3837.3; a = 13:52:17.51, <5 = -38:37:16.82). It is 
a well detached system with a circular orbit, a short orbital pe¬ 
riod (~2.5 days), and observed parameters consistent with the pre- 
main-sequence evolutionary stage. We also announce the discovery 
that it is a triple hierarchical system. 

This system was first studied in 1954 in the Cape Photo¬ 
graphic Catalog 1950 djackson & Stovl[l954b . but only the mag- 
nitude s and epoch . Then in 1984 the spectral type of F5V was de¬ 
fined dHoukl ll 982h. The syst e m is d escribed as an Eclipsing Al¬ 
gol (E A) by l Otero & Claus I d2()04l) and has a V magnitude of 
8.551 dAnderson & ErancisI 2012h . As a bright star, it was a sub- 
ject of spectroscopic s t udies of the Geneva- Copenhagen survey 
dNordstrom et al. ||2004 iHolmberg et ^l2009h . where was treated 
as a single star, however due to large brightness ratio in the visual 
we can consider some of their results reliable. The parallax and 
proper motion are also kno wn: tt = 8.43(94) ma s, = -72.49(87), 
= -44.20(78) mas yr“^ dvan Leeuwenll200% . 

In Section 2 we describe the observations used to identify 
VI200 Cen as an eclipsing and spectroscopic binary with a cir- 
cumbinary companion, in Section 3 we determine the orbital pa¬ 
rameters of the system and physical properties of the component 
stars of and compare them with theoretical isochrones. Finally, in 
Section 4 we discuss the future observations needed to fully analyse 
the system. 


2 OBSERVATIONS 
2.1 Echelle spectroscopy 

Observations with PUCHEROS instrument were carried out be¬ 
tween January and May 2012 at the 50-cm telescope of the Obser¬ 
vatory UC Santa Martina located near Santiago, Chile. The tele¬ 
scope is the European Southern Observatory (ESO) instrument, 
formerly located at La Silla. PUCHEROS is the Pontificia Uni- 
versidad Catolica High Echelle Resolution Optical Spectrograph 
developed at the Center of Astro-Engeneering of Pontificia Uni - 
versidad Catolica de Chile jlnfante et al. ll20lH : IVanzi et al.ll201^ . 
The spectrograph is based on a classic echelle design, fed by the 
telescope through an optical fibre and it covers the visible range 
from 390 to 730 nm in one shot reaching a spectral resolution of 
20,000. It is located in a stable position, in a room adjacent to the 
telescope. The science image was typically created from four sepa¬ 
rate observations, 15 minutes each, combined later into one frame. 
For the wavelength calibration we used exposures of ThAr lamps 
taken before and after the science sequence. The data reduction of 
PUCHEROS spectra was performed with the task echelle in the 
IRAE package, following standard steps: CCD reduction, spectrum 
extraction and wavelength calibration. 

The data from PUCHEROS were supplemented with spec¬ 
tra observed with the CORALIE spectrograph at the 1.2 m Euler 


telescope. CORALIE is a fibre-fed cross-dispersed echelle spectro¬ 
graph, it covers the spectral range from 380 to 690 nm reaching a 
spectral resolution of 50,000. Observations were made in a simulta¬ 
neous wavelength calibration mode, where the light of the object is 
collected by one fibre, and of the ThAr lamp by the other. The spec¬ 
tra were taken between June 2012 and July 2013. CORALIE data 
were reduced with a dedicated. Python-based pipeline djordan et al.l 
|2014 . 


2.2 Photometry 


Photometric data were obtained from the All-Sky Au tomated Sur¬ 
vey (ASAS) ( |Poimanskil2002l : |Paczvnski et al. |200^ and from the 


SuperWASP transiting planet survey ( Pollacco et al.ll200d) . ASAS 
has produced an extensive catalogue of variable stars (ACVS) of 
the southern hemispherJE In this work, we use 495 data points 
from the third stage of the survey, obtained in the V filter between 
2000 and 2009. SuperWASP is a wide-field photometric variabil¬ 
ity survey designed to detect transiting gas-giant planets around 
bright main sequence stars. The survey cameras observe bright stars 
(V~ 9-13) at high precision (1%) using a broad V + R band filter. 
3234 data points have been extracted from the SuperWASP public 
archivfl 


3 ANALYSIS 

3.1 Radial velocities and orbital solution 

Radial velocities (RVs) were measured with an implementa- 
tion of the two-dimensi onal cross-correlation technique TODCOR 
dZucker & Mazehlll994h with synthetic spectra used as templates. 
The formal RV measurement errors were computed from the 
bootstrap analysis of TODCOR maps created by adding randomly 
selected single-order TODCOR maps. The peaks of the cross¬ 
correlation function (CCF) coming from both components were 
significantly different in height, due to the large brightness ratio of 
the two stars. The PUCHEROS spectra usually had lower signal- 
to-noise (S/N) ratio (5-30) and the CCF peak from the secondary 
was not always recognized. The S/N of the CORALIE data was 
normally higher (25-60) and the secondary’s CCE peak was more 
prominent, still very low though. With TODCOR we also tried to 
estimate the intensity ratio of the two components, but due to low 
S/N of the cool secondary the results were very uncertain and we 
find them unreliable. 

The orbital solution was done simultaneously on all RV 
measurements with our simple procedure, which fits a double- 
Keplerian orbit using the Lev enberg-Marquartd algo rithm (for a 
more detailed description, see iHelminiak et aHl2009h . and allows 
for Monte-Carlo and bootstrap analysis to obtain reliable estima¬ 
tions of uncertainties. We used the photometric data first to find the 
orbital period, and check for eventual eccentricity (see next Sec¬ 
tion). The eccentricity was found to be indistinguishable from zero 
so it was held fixed to 0.0 in the further steps. In the orbital analysis 
also the period and To were held fixed. In the final orbital analy¬ 
sis we have also kept fixed the period and To. We the used values 


^ http://www.astrouw.edu.pl/asas/?page=acvs 
^ http://exoplanetarchive.ipac.caltech.edu/ 
applications/ExoTables/search.htmi?dataset 
=superwasptimeseries 
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Table 1. Radial velocity (RV) measurements of V1200 Cen, together with the final measurement en'ors (cr) and residuals from the final three-body fit (O — C). 
All values in km s~^. In the last column, 5/P denotes OUC-50cm/PUCHEROS and E/C Euler 1.2m/CORALIE observations. 


JD-2450000 

Vi 

0-1 

O-Ci 

V2 

0-2 

0 

1 

Tel./Sp. 

5714.615861 

45.958 

0.646 

-0.340 

— 

— 

— 

5/P 

5736.539995 

64.395 

0.494 

0.222 

-127.519 

2.640 

-0.757 

5/P 

5737.639889 

-67.029 

0.523 

0.394 

88.377 

4.198 

-1.718 

5/P 

5750.604835 

-62.826 

2.446 

1.482 

— 

— 

— 

5/P 

5751.584224 

74.651 

0.536 

-0.788 

-126.370 

4.324 

3.856 

5/P 

6066.642808 

47.460 

1.498 

-0.858 

— 

— 

— 

5/P 

6066.665643 

51.655 

0.841 

0.658 

— 

— 

— 

5/P 

6078.565477 

-36.335 

2.129 

1.846 

— 

— 

— 

5/P 

6080.625298 

-89.867 

0.163 

-0.009 

112.325 

1.503 

0.657 

E/C 

6081.564728 

52.113 

0.228 

-0.021 

-116.745 

1.075 

-0.456 

E/C 

6179.474281 

-26.024 

0.167 

0.035 

80.336 

0.885 

-1.105 

E/C 

6346.690592 

-12.831 

0.169 

-0.321 

67.855 

0.876 

0.573 

E/C 

6348.857536 

-55.020 

0.165 

0.255 

136.192 

1.064 

1.044 

E/C 

6349.894755 

94.865 

0.194 

0.023 

-107.687 

1.017 

-0.499 

E/C 

6397.520928 

38.353 

0.112 

-0.034 

-71.655 

0.772 

1.072 

E/C 

6398.517694 

-77.575 

0.116 

0.017 

112.000 

0.951 

-0.331 

E/C 

6497.610599 

-67.667 

0.157 

-0.152 

133.439 

0.797 

0.229 

E/C 

6498.610654 

64.361 

0.113 

0.082 

-78.099 

0.942 

0.453 

E/C 


found by JKTEBOP (see Table 3), i.e: 

Mini = 2451883.8813(19) + £;2.4828752(22) (1) 

The weight of the points was scaled according to the formal errors 
found by TODCOR. The procedure also allows for fitting for the 
differences of the zero points of different instruments, separately 
for each component. Those shifts were found to be of a similar 
value than their uncertainties, and different for each star. This can 
be explained by a mismatch between the spectrum and the template 
used, and large rotational velocities of the two stars. 

We started with a purely double-Keplerian solution with no 
perturbations, but the solutions we were getting were not satisfac¬ 
tory, with high reduced — 1600, rms of the fit about 10 km/s 
for both components, and residuals of both components correlated, 
i.e. differing from the model by a similar value. We then used a 
modified version of our procedure to look for a third, circumbinary 
body in the system. We fitted for the outer orbit’s period P 3 , ampli¬ 
tude of the inner pair velocity variations K 12 , and the base epoch 
T 3 , eccentricity es and the argument of pericentre oj. We ran the 
fitting procedure again on the whole data set, and found a satisfac¬ 
tory solution, characterised by a much lower rms - 0.98 and 1.50 
km/s for the primary and secondary, respectively and much lower 
reduced ~ 1.29. All of our spectra have barycentric correction, 
we used IRAF’s bcvcor for the PUCHEROS data and in the case 
of CORALIE data the correction is implemented in the reduction 
pipeline. Therefore, the scale of the perturbation and the outer or¬ 
bit’s eccentricity can’t be explained by the improper correction for 
the barycentre. Hereafter, following the usual convention, we will 
refer to the eclipsing pair as AB, and to the third body as C. 

We present all the measurements, together with the errors and 
residuals, in Table [T] In Table we present our results of the full 
orbital analysis. The parameters of the circumbinary orbit should 
be treated as preliminary, but the binary pair’s orbital elements are 
well constrained. The RV measurement errors were initially scaled 
in such way that the final reduced was close to 1 , and the fit it¬ 
self was not affected (weights not changed). The errors were found 
by a bootstrap analysis (10000 iterations). In such a way we take 
care of the possible systematics and obtain reliable uncertainties of 


Table 2. Results of the RV analysis and the orbital parameters of 
VI200 Cen. 


Pai'ameter 

Value 

± 

AB eclipsing pair 


Ki [km s“^] 

78.23 

0.37 

K2 [km s“l] 

126.0 

1.1 

V'y [km s~^] 

10.92 

0.94 

ai2 sini [Rq] 

10.026 

0.058 

ei2 

0.0 

(fixed) 

<? 

0.6208 

0.0062 

Ml sin^ i [M0] 

1.352 

0.027 

M2 sin^ i [M0] 

0.839 

0.012 

Outer orbit 


P 3 [d] 

351.5 

3.4 

Ts [JD-2540000] 

5358.2 

8.0 

Ki 2 [km s“i] 

18.37 

0.11 

63 

0.42 

0.09 

OJs 

156 

7 

as sin is [AU] 

0.538 

0.033 

/[Mq] 

0.099 

0.025 

Ms{is = 90°) [M©] 

0.662 

0.066 

Other fit parameters 


E/C-5/Pi“ [kms-l] 

0.98 

0.31 

E/C-5/P2“ [kms-l] 

5.32 

2.00 

DoF 


20 

rmsi [km s“^] 

0.68 

rmsi [km s“^] 

1.40 


“ E/C—5/P is the difference in spectrograph zero 
points measured for each component separately. 


the final parameters. Our model, separated into the Keplerian and 
perturbation components, is presented in Figure[T] 

It is worth noting that the mass function /(M 3 ) was calculated 
from the full formula: 


/(M3) = __ 

(Mi+M2+M3)2 


( 2 ) 


and is quite large. The minimum mass of the third body (for is = 
90°) is a substantial fraction of the total mass of the eclipsing pair. 
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_ I _I_ I _I_ I _I_^_ 

5800 6000 6200 6400 

JD-2450000 


Figure 1. Three-body orbital model of V1200 Cen based on the RV mea¬ 
surements from PUCHEROS (circles) and CORALIE (triangles). Filled 
symbols are for the primary component, and open ones for the secondary. 
Black dotted line on panel a) marks the systemic velocity Vj. a) Keplerian 
orbit of the AB pair as a function of the orbital phase, with the perturbation 
removed; b) residuals of the full model as a function of phase; c) pertur¬ 
bation from the third body as a function of time, Keplerian orbit removed; 
d) residuals of the full model as a function of time. 


so the popular approximation (Mi + M 2 + Mz)^ — (Mi + M 2 )^ 
is not valid. We checked for tertiary eclipses in the residuals of 
the light curve models (see next Section) and we found no obvious 
evidence for them, which means that the outer orbit’s inclination 
is different from 90°. We have also failed to identify a tertiary 
peak in the CCF, which means that the tertiary component is sig¬ 
nificantly dimmer than the secondary. All this is consistent with 
for example a 0.7 M© star on a 70 — 75° orbit. It is also not ex¬ 
cluded that the component C may itself be a binary composed of 
two similar, lower-mass stars, whic h together exceed the mass of 
B. This would be in agreement with iTokovininI ( l2014h . who found 
no hierarchical triples with outer orbits shorter than ~1000 days. 
The mechanism leading to a formation of such double-double sys¬ 
tem, with short inner orbital periods and outer period shorter than 
1000~ d, assumes non -aligned outer and inner orbits jWhitworthl 
1 200 ll ; Irtakovininl [ 20 1 4) . It is thus possible that in V1200 Cen the 
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3 

3 
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- 0.05 
I 0 
° 0.05 
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.S 8.6 

3 

-hJ 

'S 
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8.8 


- 0.05 
I 0 
° 0.05 



- 0.2 0 0.2 0.4 0.6 0.8 

Phase 


Figure 2. Observed ASAS (top) and SuperWASP (bottom) light curves with 
JKTEBOP models over-plotted. Lower panels show the residuals of the fit. 


mutual inclination between the outer and inner orbit is large, mak¬ 
ing the Kozai mechanism possible (Kozai 1962). 


3.2 Light curve modelling 


Both ASAS and SuperWASP light curves were fitted using the JK¬ 
TEBOP code (Southworth et al. 2004a,b), which is based on th e 
Eclipsing Bi naries Orbit Program (EBOP; [Popper & Etzellll9^ . 
and P HOEBE j Prsa & Zwittej200^ - an implementation of the WD 
code jWilson & Devinn^ 197lh . JKTEBOP determines the optimal 
model light curve that matches the observed photometry and re¬ 
ports the parameters obtained for the model. The parameters de¬ 
rived include the period P, time of minimum light To, surface 
brightness ratio, relative sum of the radii (i?i + R 2 )/a, ratio of 
the radii R 2 /RI, inclination i, eccentricity e, and argument of pe- 
riastron uj. The routine also takes into account the effects of limb 
darkening, gravity darkening and refle ction effects. We ado pted the 
logarithmic limb darkening law from IClaret et"^ l l2000l) . To ob¬ 
tain the limb darkening coefficients we considered the ASAS light 
curve to be in Johnson’s V passba nd. For the SuperWAS P filter we 
interpolated the coefficients from IPollacco et al.l l l2006l) . We also 
checked the ASAS data for 0-C timing variations, and we did not 
find any significant ones. Thus we did not correct the light curves 
for the light time effect induced by the third body. We used PHOEBE 
to obtain the temperature of the secondary star, u sing the temper¬ 
ature of the primary star known previously from IHolmberg et al.l 
(l2009h . We did not use the latest Geneva-Copenhagen survey re¬ 
sults, as they rely on infrared data, which in the case of VI200 Cen 
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Table 3. Pai'ameters obtained from the the ASAS and SuperWASP (SW) LC analysis. The adopted values are weighted averages. 


Parameter 

ASAS Value 

± 

SW Value 

± 

Adopted Value 

± 

p[d] 

2.4828778 

0.0000043 

2.4828752 

0.0000025 

2.4828752 

0.0000022 

To [JD] 

1883.8789 

0.0031 

1883.8827 

0.0024 

1883.8813 

0.0019 

in 

81.9 

+2.8 

-1.3 

81.6 

+ 1.6 
-1.3 

81.8 

+ 1.4 
-1.2 

ri 

0.137 

+0.014 

-0.015 

0.138 

+0.025 

-0.034 

0.137 

+0.014 

-0.015 

r2 

0.107 

+0.024 

-0.039 

0.110 

+0.038 

-0.026 

0.109 

+0.022 

-0.025 



Figure 3. Results of the eiTor analysis performed with JKTEBOP with Monte-Carlo on the ASAS (red) and residual shifts on the SuperWASP data (green). 
Plots present the distribution of consecutive solutions on the ri + r 2 vs. i (left) and k = vs. i (right) panels. Black stars with error bars correspond to 

the adopted values with their la uncertainties. 


can be affected by the secondary and tertiary star. The value of 
the mass ratio was found in the previous orbital analysis. Figure|2] 
shows the observed ASAS and SWASP light curves together with 
their models. 

Reliable uncertainties in case of ASAS data were calculated 
with the Monte-Ca rlo method (10000 runs) and with residual-shifts 
(ISouthworthl2008l) in case of SWASP data. Figurej^shows the dis¬ 
tribution of ri + r 2 and k as the function of the inclination angle. 
The usual correlation is clearly visible. The results of the LC analy¬ 
sis with JKTEBOP are presented in Table[^ We found that the ratios 
of radii (k) and effective temperatures are the most uncertain val¬ 
ues, vastly contributing to errors of such physical parameters as 
absolute radii or luminosities. This is due to large scatter of the 
photometric data and low contribution of the star B (12% in V). 
One way to constraint k would be to use intensity ratios from TOD- 
COR, but as we mentioned the S/N of the spectra was too low. Data 
of a much higher S/N are needed, optimally taken in IR (both pho¬ 
tometry and spectroscopy), where the secondary’s contribution is 
larger. 


3.3 Kinematics 

Using the known parallax and proper motion l lvan Leeuwen 11200% 
and our value of the systemic velocity n-y, we have calculated the 
galactic velocities: U = —36.7 ± 3.3, V = —21.8 ± 3.6 and 
W = —1.8 ± 0.6 km s“^ (no correction for the solar movement 
has been done). These values put V1200 Cen in the Flyades moving 


group JSeabroke & Gilmorel l2007l : lZhao_et^ 200^, which sug¬ 
gests the age of -^625 Myr. Fiowever. iFamaev et al l l l2008l) have 
shown that about half of the stars that reside in the same area in the 
velocity space as the Hyades group, actually does not belong to it. 
Nevertheless, V1200 Cen seems to be a young system belonging to 
the thin galactic disk. 


3.4 Absolute parameters 

The absolute dimensions and distance were calculated using JK- 
TABSDIM with the results obtained from the radial velocity and 
light curve analysis. The parameters with their respective uncertain¬ 
ties used for the input were the velocity semi amplitudes (km s~^), 
period (days), orbital inclination (degrees), fractional stellar radii 
(i.e. in units of the orbital major semi-axis), effective temperatures, 
and apparent magnitudes in the B and V filters. We did not use the 
available JHK photometry as it could have been affected by the 
third star. For the interstellar reddening we considered E(B-V)=0, 
and we observed no major differences for the distances obtained 
in the two filters, both in agr eement with the valu e of the parallax 
8.43±0.94 mas (119±13 pc: Ivan Leeuwen ll2007h . We present the 
the physical parameters of V1200 Cen in Tablej?) 


3.5 Evolutionary status 

In Figure|4]we present our resul ts from Table l^plotted over theoret¬ 
ical stellar evolution models of isies 3 1200^ and from Yonsei-Yale 


© 2014 RAS, MNRAS 000,[T]{8] 


























6 J. Coronado et al. 


Table 4. Physical parameters of V1200 Cen obtained with JKTABSDIM on 
the basis of spectroscopic values of and [Fe/H]. 


Parameter 

Value 

± 

a (R©] 

10.13 

+0.07 

-0.06 

Ml [Mq] 

1.394 

0.030 

M2 IMq] 

0.866 

0.015 

Ri [R©] 

1.39 

+0.14 

-0.15 

R2 IRq] 

1.10 

+0.22 

-0.25 

log gi 

4.30 

-0.09 

+S4S 

log 92 

4.29 

-0.18 

+0.20 

[Fe/H] 

-0.18“ 


Teff,l IK] 

6266“ 

94 

Teff,2 IK] 

4650*’ 

900 

Ll IlogjL/L©)] 

0.42 

+0.09 

-0.10 

L 2 Ilog(L/L©)] 

-0.29 

+0.38 

-0.40 

d [pc] 

98 

11 


“ From iHolmberg et ^ i2009l) . 

^ From temperature ratio obtained with PHOEBE. 

rY^: lYi et alfcOOlh set for ages of 6, 10, 15, 20 and 625 Myr. For 
the Y^ we also include the 3 Gyr isochrone. 

The values from Table|4]are marked as black points. The two 
most precise values we have are the mass and temperature of the 
primary. One can see that it is too cool for a main sequence ob¬ 
ject, so it is either evolved, as the 3 Gyr isochrone suggests, or at 
its pre-main-sequence (PMS) stage. However, both Y^ and Siess’ 
models predict a drastic change in temperature and radius between 
10 and 20 Myr. The primary’s radius is much too small for an 
evolved star and agrees with the late PMS, meanwhile the tempera¬ 
ture clearly suggests younger age, inconsistent with the radius. We 
also note that secondary’s parameters are much better reproduced 
by PMS models, especially its radius, which is much larger than 
main-sequence objects of the same mass, as expected for stars that 
are still evolving onto the main sequence. 

The comparison of our results with the model predictions in¬ 
dicates that V1200 Cen is a pre-main-sequence system, but no 
age is fully consistent with the data, and the resulting distance is 
only in a fair agreement with results from the Hipparcos. Un¬ 
fortunately, the uncertainties are very large, especially for the sec¬ 
ondary component, so our conclusions cannot be treated as final. 
The strongest constrain we have comes from the primary’s tem¬ 
peratu re, which was derived spectroscopically by iHolmberg et al.l 
( 1200911 . We find almost exactly the same temperature - 6263 K - 
using the observed B — V c olour (0.475 mag) and calibrations by 
ISekiguchi & Fukugital ( 120000 . One has to remember that the Holm- 
berg’s temperature was derived under the assumption that the star 
is single. Thus, it is still possible that the true Tsff is higher. We 
run a series of tests to find the temperatures that give the best agree¬ 
ment with the Hipparcos parallax, assuming their ratio to be the 
same as found by us with PHOEBE. We found the best match to 
the observed distance for much higher temperatures of 6900 and 
5120 K for the primary and secondary respectively. The resulting 
effective temperatures and related luminosities are summarised in 
Table|5] We plot them on the Figure|4]with red symbols. 

One can see that the new higher values of radiative parameters 
fit the main sequence or even 20 Myr models, and both radii still 
agree with the 30 Myr isochrone. The age of 30 Myr also happens 
to be the time scale of circularisation of the system’s orbit. All in 
all, we get a self-consistent model of a 30-625 Myr old, young mul¬ 
tiple. As an attempt to distinguish between high and low tempera- 


Table 5. Radiative parameters of V1200 Cen obtained with JKTABSDIM by 
fitting the temperature scale to match the Hipparcos distance. 


Parameter 

Value 

± 

Teff.l [K] 

6900 

100“ 

Tea,2 IK] 

5120*’ 

900“ 

Ll [log(L/L©)] 

0.59 

+0.09 

-0.10 

L2 [logjL/L©)] 

-0.13 

+0.36 

-0.37 


“ Uncertainty assumed. 

^ From temperature ratio obtained with PHOEBE. 

ture scale we also run a simplified s pectral analysis with the S pec- 
troscopy Made Easy package (SMEt lValenti & Piskunovll996l) . We 
took the CORALIE spectra, shifted by the measured velocity of the 
primary and stacked them together. We run the SME on a portion of 
spectrum spanning 6190-6260A, and keeping the [Fe/H], logjg) 
and Vrot on values from or expected from Tab.|4] The secondary’s 
contribution was treated as a contribution to the continuum, con¬ 
stant across the wavelength range. We obtained Te//,i ~ 6000 K, 
favouring the lower temperature scale, and younger ages (Fig. 0. 
We are aware of the fact that such analysis is affected by the sec¬ 
ondary, and we do not treat it as a proof for the correctness of 
Holmberg’s temperature, but we find it unlikely to be off by almost 
1000 K, between spectral types FO and F8. 


4 CONCLUSIONS 

We report the discovery that the eclipsing binary VI200 Centauri 
is a triple, likely a member of the Hyades moving group, but the 
largely inflated secondary’s radius may suggest that the system may 
be PMS around 30 Myr. There are few such objects known so far, 
however they are very important for calibrating stellar evolution 
models at young ages, where stars are changing rapidly as they 
evolve onto the main sequence. Analysis of ASAS and SuperWASP 
light curves combined with radial velocity measurements allowed 
to obtain absolute parameters of the system. Through our analysis 
of radial velocities and orbital solution we determined the presence 
of a third companion, in a wide orbit. We reached a good precision 
in mass determination (2.2 and 1.7%), but other parameters (radii, 
temperatures, luminosities) are not that well established. Further 
analysis of the system allowed us to compare our results with stellar 
evolution models, obtaining an approximate age of 30 Myr. Despite 
its brightness, data of higher S/N are required to better constrain 
physical parameters of the system, especially temperatures and ra¬ 
tio of the radii. Possible solution would be to obtain photometry 
and spectroscopy in the IR, where components B and C contribute 
more than in visual. The secondary eclipse would be deeper, and 
the influence of the third light should be detectable. Also, if RVs of 
the star C were measured, a full dynamical solution of the system 
could be obtained, including mass of the third star and the inclina¬ 
tion of its orbit. This would also allow for comparing the isochrones 
with three stars, thus constrain the age and evolutionary status even 
better. With its probable age, VI200 Cen is an important object to 
study the tidal and third-body interactions young binaries. 
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V1200 Cen: a bright triple in the Hyades moving group 1 








Figure 4. Mass vs. radius, log(L) and log(Tp fF). The left co lumn shows our results and the isochrones from for ages of 6, 10, 15, 20, 625 Myr and 3 Gyr. 
The right column shows the same but using the lSiesj i2Q0Ql) isochrones (without 3 Gyr). Black symbols are for our results from Table[4]and red symbols for 
the results from TableO 
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